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I. INTRODUCTION 

Analytic forms for flavor transition probabilities of 
neutrino oscillations in matter can facilitate studies of 
measurement sensitivity afforded by proposed new ex- 
perimental facilities [IH5] . Of particular current interest 
is the transition probability for subdominant v e appear- 
ance, — > v e , using beams propagating over long 
baselines through terrestrial matter. Indeed, the size of 
the neutrino mixing angle #13 has been dramatically clar- 
ified during the past year as the result of a number of 
recent experimental measurements [BHlOj. 

A determination which is representative of the new 
level of precision is reported by the Daya Bay reactor ex- 
periment: sin 2 26»i 3 = 0.092 ±0.016 (stat) ±0.005 (sys). 
With #13 ~ 9°, the study of CP violation and mass hier- 
archy in the neutrino sector can proceed in earnest. An- 
alytic investigations in this new era can benefit from the 
availability of convenient analytic forms for Viy^ — > i/ e ) 
which are accurate to within a few percent for neutrino 
baselines through the Earth's mantle. This work provides 
such a probability expression for — > 24 oscillations in 
a constant-density matter field. 

A number of exact derivations for neutrino propagation 
with oscillations among three active flavors in constant- 
density matter have appeared in the literature over the 
past several decades [TTHTB] . In general these formula- 
tions are complicated and do not readily yield insights. 
A degree of clarity is achieved with the exact probability 
expressions of Kimura, Takamura, and Yokomakura in 
which matter effects are disentangled from CP violation 
effects [ini HZ] . Their formulation has been extended to 
include, for example, nonstandard interaction matter ef- 
fects [18l[l9]. Nevertheless, a desire for more transparent 
formulations has led to the development of various ap- 
proximation expansions; a review with comparisons can 
be found in Ref. [20]. A treatment which incorporates 
the magnitude of #13 as recently measured into a pertur- 
bative framework is presented in [21j . 

For the v e appearance probability Piy^ — > u e ) the 
three-term formula of Cervera et al. [55] (see also [53]) is 
frequently utilized in analytic studies. This formula has 
the form of a perturbative expansion in terms of the small 
mass hierarchy ratio a = Am^/Amjj ~ 1/32. For prop- 
agation through the Earth's crust such as occurs with 



the T2K (295 km), MINOS (735 km), and NOuA (810 
km) baselines, the formula of Ref. [55] is adequate for 
most purposes. However for terrestrial baselines which 
exceed the proposed LBNE baseline of 1300 km, such 
as the "bimagic" 2540 km baseline [2] and the "magic" 
7500 km baseline [TJ [5] which have significant pathlength 
through the Earth's mantle, more accurate formulations 
are desirable. The Cervera et al. formula can be written 
as follows: 
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In the above expression and throughout this paper, we 
use A = Am^i/ (<iE„) for the atmospheric oscillation 
phase at baseline £. The symbol A refers to the matter 
potential A = ±(2y/2Gpn e E u ) / Am^ where Gp is the 
Fermi coupling constant and n e is the electron density 
in matter. The sign of A is determined by the sign of 
Am§i and choice of neutrino or antineutrino propaga- 
tion. As a matter of convention, the formalism of this 
work refers to neutrino propagation (the + sign) and as- 
sumes the normal mass hierarchy for the neutrino mass 
eigenstates (Amf^ positive). We also use the compact 
notations Sjj = sin0y and cy = cos0,j with i,j = 1,2,3. 

We present a formulation of the v e appearance prob- 
ability V^Uy. —> v e ) which retains the convenient pertur- 
bative form of Eq. ([I]) in the leading three terms of a six- 
term expansion. Our derivation uses conventional non- 
relativistic quantum mechanics methods to construct an 
approximate form for the time evolution operator for neu- 
trino states in flavor basis. Amplitudes of useful precision 
are thereby implied for all oscillation transitions accessi- 
ble to neutrinos of three active flavors. In this work how- 
ever we focus upon — > v e oscillations. Our approach is 
amenable to augmentations such as obtaining oscillation 
amplitudes with inclusion of selected non-standard inter- 
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action (NSI) matter potentials, but such developments 
are left for a future work. 



II. OUTLINE 

The paper proceeds as follows: We define some con- 
venient notations and proceed straightaway in Sec. |III| 
to report our result, stating the three amplitudes (Sec. 
Ill A[ ) and the consequent six-term oscillation probabil- 



ity (Sec 



III B ) which comprises our formulation for 
; ) in a constant density medium. We then 



show that our leading terms have resemblances to the 
probability terms of Eq. . The remaining three terms 
of our expression entail relatively small contributions to 
Viy^ — > v e ). The extent of variations is illustrated using 
plots of our six-term probability versus Eq. for E v 
between one and ten GeV at the proposed LBNE, the 
bimagic, and the magic baselines (Sec. HI C| . 

The derivation of our V(i>n — > Ve) formula is given in 
Sec. 



IV 



In brief, the conventional three-flavor Hamil- 
tonian in flavor basis is transformed to propagation ba- 
sis, whereupon it is re-phased and then separated into 
an "unperturbed" piece plus an interaction potential of 
relatively small coupling strength (all elements propor- 
tional to sin 2#i2 • a). The latter is used to set up an 
interaction picture and the time evolution operator for 
neutrino states is constructed in a heuristic way in that 
picture. A key step is the exponentiation of the inter- 
action potential into exp(— zV/(t)); the latter can be ex- 
pressed as a matrix identity which has similar form to 
the well-known identity for exp(—iJ y 9) for j = 1 angular 
momentum. Our heuristic form of the evolution operator 
is then transformed back into propagation basis and fi- 
nally into flavor basis, yielding our approximate solution 
for lJ( a \ In that basis, its matrix elements coincide with 
the various three-flavor neutrino oscillation amplitudes. 



The v e appearance amplitude is given by U\ 



(«) 



III. FORMULAE FOR v u 



OSCILLATIONS 



A. Conventions and notations 

In expressions to follow we refer to the vacuum oscil- 
lation length £ v = E^/ Am 31 , and we write 

A = Am 2 31 £/(AE U ) = -L. and V e = . (2) 

Convenient quantities are the neutrino mass hierarchy 
ratio a ~ 1/32 and the scaled forms for the hierarchy 
parameter a' and a" defined as 



sin 2#i2 • a , and 



(3) 



a" = (l- 3c 2 2 ) ■ a . 

In our expressions, the mixing strengths involving # 13 
are often accompanied by the factor (1 — s 2 2 a), hence we 
define 

sin2# 13 = (1 — s\ 2 a) ■ sin2<9 13 , and 

cos 20i3 = (1 - s\ 2 a) ■ cos 2d 13 . (4) 



B. Oscillation amplitude of three terms 

The amplitude for v e appearance at baseline £ from an 
initial beam which propagates through matter of con- 
stant density, can be expressed as a sum of three terms: 

A(v fl ^v e ) = T 1 +T 2 +T 3 . (5) 

The individual Tj terms are the following: 

T 1 = H)sm2^ 3S23 - S ^-e-^, (6) 
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and 
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The factor N which appears in oscillation phases in Eqs. 
^ and (|8| is defined by 
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In Eqs. and we use 



The quantity Fa appearing in Eq. (|8| is 



F A 



cf 3 (l - s\ 2 a) - (cos 26»i3 - -4) 
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The variables designated in ([9]), jlOj ), and ( 11 ) are ones 
which arise naturally in the derivation of Sec. |IV| To fa- 
cilitate comparison with Eq. Q we also define as conve- 
nient quantities D and A': 



D = 4£ V N ~ \1-A\, 



A' = G£ = A [(1 + A) + a"} 
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Thus the oscillation phases N£ and rj£ can be written as 
(D/M v )£ = DA and as (a'/U v )£ = a' A respectively. 

C. Oscillation probability of six terms 

The i/fi — ¥ v e oscillation probability is constructed from 
\A{Vfj. v e )\ 2 - Referring to Eq. |5]), we express the 
result as a sum of six terms: 

P{ V(l -> v B ) = I Ti| 2 + (TiT 2 * + T*T 2 ) + 



\T 2 \ 2 + {T 1 T*+T*T 3 ) 

(t 2 t; + t;t 3 ) + \t 3 \ 2 
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We proceed to construct the individual probability terms 



which appear in Eq. (14). Squaring the amplitude of Eq. 
Q we obtain 



\T^ = (sin20 13 ) s 



2 2 sin 2 (DA) 



(15) 
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Similarly, we obtain from Eq. ([T]) 

\T 2 \ 2 = c\ 3 c\ 3 ■ sinV) = c\ 3 c\ 3 • sinVA). (16) 

Considering the T\T 2 cross terms, we write 
{T{r;+T*T 2 ) = 

sm{N£) 
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(17) 



Upon introducing the notations of (12) and (13) and re- 



ducing the bracket expression, we obtain 
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There are two more sets of cross terms to consider. We 
first determine (TiT 3 * +T X *T 3 ): 

(TxT 3 * + T*T 3 ) - - 2 sin 29 13 sin 26 13 s 2 23 F A 

a /VAN sin 2 (L>A) (19) 
• Sm (-yj-— 

The remaining cross term is (T 2 T* + T*T 3 ). Now T 2 T* 



is 
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Adding T%T 3 to T 2 T 3 and extracting the common factors, 
we obtain 



(t 2 t; + t;t 3 ) = 

i sin26»i 3 ci 3 sin26> 23 ■ sin(r^) ■ sin' . 
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Equation (21) reduces to 
{T 2 T;+T*T 3 ) = 

sin 26*i 3 C13 sin 2# 23 • sin(a'A) • sin 2 
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The final term is |T 3 | 2 . Referring to Eq. (18]) we obtain 



|T 3 | 2 = sin 2 20 13 4, -sin' 
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The sum of the six probability terms of Eqs. (15), 
(16), (18), (19), (22), and (23) comprise our rendering of 



the probability for v e appearance in an initial beam. 
For convenience we state the total probability as a single 
expression in the conclusion (Sec. [v| of this work. 



D. Comparison to the Cervera et al. probability 



(|14| for the — > v e oscillation probability, the 



In Eq, 

three leading terms are | T\ | 2 , (TiT 2 * + T{T 2 
for which explicit expressions are given by Eqs. (15), 



(18), and (16) respectively. These terms have a clear 



resemblance to the corresponding three terms of the for- 
mula given in Eq. (JlJ, however our extended perturbative 
form gives rise to certain modifications. In the two lead- 
ing terms, our variables sin26'i 3 and D replace sin26'i 3 
and (1 — A) respectively. For the second and third terms, 
the correspondence between our result versus the Cervera 
et al. expression can be seen by invoking small angle ap- 
proximations and by recognizing that A' ~ A and that 
1.0. 



°13 

Equation ( 14 ) contains three additional terms which 
arise from the presence of the very small T 3 ampli- 
tude. These terms go beyond the level of accuracy in- 
tended with Eq. Q. We find that these extra terms 
contribute amounts to the probability of less than one 
percent for terrestrial baselines accessible to accelerator- 
based oscillation experiments. Thus there is justification 
for neglecting these relatively complicated higher-order 
terms. The improved accuracy afforded by our formula- 
tion arises in the main with the refinements introduced 
into the three probability terms already present in the 
Cervera et al. formula |22j . rather than in the extra 
terms. 

To illustrate the level of improvement, we show the 
— > v e probability for neutrinos of energies between 1.0 
and 10 GeV, propagating through the Earth for three 
baselines of interest to future experimentation. The 
nominal value reported by the Daya Bay experiment, 
sin 2 26*13 = 0.092 [9], is used throughout, and the normal 
mass hierarchy is everywhere assumed. Figure [T^, shows 
the v e appearance probability for the Fermilab to Home- 
stake baseline of 1300 km as envisaged for the Long Base- 
line Neutrino Experiment (LBNE). At this baseline neu- 
trino propagation is entirely through the terrestrial crust, 
and a uniform density of 2.72 g/cm 3 is assumed. Our 
six-term formula (solid curve) agrees with the more ap- 
proximate probability (dashed curve) fairly well, however 
a small reduction in v e appearance is indicated through- 
out the peak oscillation region. The difference between 
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the predictions becomes negligible at shorter terrestrial 
baselines. 
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FIG. 1. Probability for v e appearance from accelerator 
neutrino beams propagating in constant density matter, for 
the (a) LBNE and (b) bimagic baselines. In each plot our 
six-term transition probability (solid curve) is compared to 
the three-term probability of (IT]) (dashed curve). Our for- 
mula predicts the probability to be somewhat lower through- 
out the E v region of peak oscillation. This disparity is more 
pronounced at the longer baseline. 



For baselines longer than LBNE the reduction in v e 
appearance probability becomes more pronounced. Fig- 
ure [TJa compares expectations at the bimagic baseline of 
2540 km. At this baseline, there is propagation through 
the Earth's mantle as well as the crust, giving rise to a 
mean density of 3.2 g/cm 3 . In the vicinity of the oscilla- 
tion peak at ~ 4 GeV our result falls below the Cervera 
et al. probability by several percent; this trend that per- 
sists at higher energies. The same general trend has been 
shown by other improved approximation forms - see for 



example, Fig. 3 of Ref. [H]. 

With even longer baselines, the MSW resonance in the 
mantle |24j greatly enhances the — > v e oscillation. 
Figure [2] compares Eqs. (14) and ([T]) at the so-called 
magic baseline which occurs for propagation paths in the 
neighborhood of 7500 km. Here the propagation is pre- 
dominantly mantle traversal and the mean density is 4.3 
g/cm 3 . The resonance-driven appearance probability is 
nearly 50% at its peak, consequently Fig. [2] is plotted 
with a distinctly larger abscissa range than is used in 
Figs. [lji,b. Our probability formula exhibits the same 
shape as predicted by Eq. ([!]) over most of the E v range. 
However it shows the appearance probability to be over- 
estimated by Cervera et al. throughout the region of the 
main oscillation peak. Figures [JJ and [2] indicate the ex- 
tent to which the six-term probability of this work may 
offer improved accuracy for long baseline oscillations. 



1.0 



><0.8 



I = 7500 km 
S C P -0 



■ This work 
Cervera ef al. 




4 6 8 

Neutrino Energy [GeV] 



10 



FIG. 2. Probability for — > v c oscillation at the magic 7500 
km baseline through the terrestrial mantle. Expectations for 
our six-term formula (solid) and for the Cervera et al. result 
(dashed) are shown; a mean density of 4.3 gm/cm 3 is used for 
both. The effect of the MSW resonance in the Earth's mantle 
is apparent in both curves; our formula gives a probability 
which is lower by ~ 8% in the vicinity of the first oscillation 
peak. 



Figurej3] compares |T\ | 2 of Eq. ( 15 ), with the lead term 
of Eq. (nT. For the Cervera et al. result , the leading 
term gives the entire probability contribution at 7500 km 
whereas the second and third terms, which carry the CP- 
violating phase and parameters from solar-scale mixing, 
have negligible probability at the magic baseline. 

Figure [3] shows our |Ti| 2 term corresponds to a proba- 
bility which is nearly 20% below the Cervera et al. pre- 
diction in the maximum oscillation region. This dis- 
agreement is partially alleviated by modest contribu- 
tions, mostly of positive sign, which in our formula arise 
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from the terms (Ti^ + T^) and |T 2 | 2 . The three other 
terms in our formula, namely those which involve the T 3 
amplitude, contribute an amount which is only ~ —0.002 
throughout the interval 2.0 < E v < 10.0 GeV. 
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FIG. 3. Comparison of the leadi ng t erm in Viy^ — > v e ) 
from the formula of this work (Eq. ( |15[ )) to the correspond- 
ing term in Eq. |l]), for the 7500 km baseline, the leading 
terms exhibit similar shape dependence with respect to neu- 
trino energy, however |Ti| 2 of our formulation gives a lower 
probability through the region of the MSW resonance. 



IV. DERIVATION OF A(y^ -> v e ) 

A. Hamiltonian in flavor basis 

For neutrino propagation in vacuum, the Hamiltonian 
in the basis of three mass eigenstates Vi (i =1, 2, 3) is 

Hf =diag(0,^L ^} = £ -diag(0,a,l). The 

transformation from mass basis to neutrino flavor- 

basis {if — e, /i, t) is provided by the unitary 

mixing matrix 



(24) 



The Dirac-type CP-violating phase, Sep, can be conve- 
niently incorporated into Uf m i X ) by including auxiliary 
matrices is CP within the standard factored form: 

U (mix) = R^Baa) ■ h CP • 4(0X3) ■ t-a OP • Mfa) (25) 



where Ig CP = diag(l, 1, e i5cp ) and l\ 



diag(l,l,e lAcp ). In Eq. (251 the atmospheric and so- 
lar mixings are accounted for via the rotation matrices 



OCP 



■fix (^23) anc! £3(0x2), and the product ls CP -R 2 {9i 3 )-l-s ct 
carries the CP-violating phase: 



■ Rn ■ 



-Set 



c 13 si 3 e- l5cf 
1 
-s 13 e l5cp C13 



The effective wave equation for vacuum propagation of 
flavor states is 



(26) 



Here, the vacuum Hamiltonian in flavor basis is given by 
the unitary transform of , 



-^0 = yRi^s C p^2^—Scp^3j H Q yR 3 Is ap R 2 I-s ap R 1 

(27) 

Since l_{ OJ> (Ij OP ) commutes with R 3 (R^), and since 



_S op Hq % cp = Hq ', Eq. (|27|) can be simplified: 
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c 12 a ~ 2 Sl3Q; 
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Upon inclusion of the MSW matter interaction 



H, 



:<p) 



diag (14,0,0) , 



the total Hamiltonian including matter effects as well as 
CP violation can be written in flavor basis as 



Hto = ( fii 



(h {23) +H iv) 

1 T -"matter 



) (is CP Rl) ■ (29) 



Thus the full wave equation in flavor basis is Eq. (26) 
with defined by Eq. (29) replacing on the 

right-hand side. 



B. Hamiltonian in propagation basis 

Neutrino propagation is usefully re-cast by transform- 
ing to the propagation basis. The latter basis is defined 
via 

= l_ Scp R T l {9 23 )v^ , ^=R x t 5ap iP>\ (30) 

Multiplication of the wave equation in flavor basis from 
the left by I_s ap R[ yields 



i-z>)(t)=#Wi?(p) . 
at 



(31) 
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where = ( H, 



(23) 



flatter ) i s the effective Hamilto- 



nian in the propagation ba sis. The matrix is nearly 



identical to H, 



(23) 



of Eq. ( 29 1 but includes the matter 

(23) x 



term V e — Aj2l v added to the element (Hq )u- Hamil- 
tonian is a real- valued, symmetric matrix devoid of 
the CP- violating phase 6c p- Of course the CP-phase 
reappears when one transforms from propagation basis 
back into flavor basis. 

The matrix can be "re-phased" . That is, we per- 
form an algebraic manipulation leading to removal of a 
term proportional to I, which merely contributes an over- 
all phase to the oscillation amplitudes. Specifically we 
subtract (and also add, but then discard) the following 
diagonal matrices to H^ p >: 



r 2 - 


1 


1 


t 2£y 


2 










s{ 2 a 


fc\ 3 






V 2 " 


-T 



(A - cos26>i 3 )I 
1 



and 
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At 



-s\ 2 a. 



Upon extraction of a factor ~ we obtain 
1 



-(cos 20i3 — A) c\ 3 a' sin2#i 3 

c 13 a' - [(1 + A) + a"} -s 13 a' 

sin 26>i 3 -Si3a' +(cos 26 13 - A) 



= —G. To represent the other elements of in 



Using the variable G defined in Eq. (j 1 0|) , we identify 

y(p) 

! 22 

a compact form, we define: 

1 , 1 



Q= —(cos 20i3 -A), /=_ s in20 13 , 



(33) 



The full Hamiltonian in propagation basis can then be 
written as 

-Q a f 
= I a -G b ) . (34) 
/ b +( 



C. Formulation in an interaction picture 

We separate into an "unperturbed" part, Hq P \ 

plus an interaction potential, V comprised of elements 
proportional to a': 



= H, 
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(35) 



Wave equation (31 ) can then be re-cast into an Interac- 
tion picture: 

VW(t) = e^ P Vp)(i), 0ip)(t) = e -4 Wt #)(f) , (36) 
so that 

i^pW(t) = V I -0M(t) (37) 

where 

V I {t) = e ik ° h -Y-e- ift « h . (38) 



Our approach to Eq. (37) is to construct a heuristic 
functional form JJ^ estS> (t, 0) which serves as an estimator 
of the time evolution operator in the Interaction picture: 

^(t)-; 0)-z7«(0). (39) 

The wave equation which defines our f/j es *^ (t i 0) is 

T (est) u n ^ z^q-j 



i~U\ est \t,Q)=VAt)-U\ esz \t,Q). 
at 



where t is variable but V/(t) is set to its value at the final 
baseline distance t = I. 

To obtain V/(i) we require the matrix representation 
(in propagation basis) of the unitary operator forms 

exp(±iHQ P h), where H^f 1 has the elements given in Eq. 



( 35 1 . Considering the series expansion 



W = e ik ° h = 

n=0 



(41) 



it is readily seen that neither the middle row nor middle 
column of matrix mixes with other elements. Then 
the matrix W has the form 








w 13 





e -iGt 





W31 





^33 



(42) 



Thus we may work with the reduced 2x2 matrix 

- f& x - Qo z (43) 



Mp) _ 

n R — 



-Q I 
f +Q 



where a x<z are the Pauli spinor matrices. We write 
H ip] = N ■ a , where N = (/,0,-Q), with \N\ = 
y/f 2 + Q 2 (see Eq. (9)). The unit vector n = N/N = 

l2 \-l/2 



(/, 0, —Q) serves as the axis-of-rotation in 

(p) 



(f + Q 2 

the reduced spinor space. With Hr' — Nn ■ <?, and reo 

(44) 



ognizing that t — I in natural units, we write 



e iH^'(t=l) _ e ih-3{Nl) _ e ih-34 
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where <fi = N£ designates the rotation angle. With n 
(n x , 0, n z ), the spinor identity is 



2™'^ =Icos(j) + ia- 



cos <p + m z sin cp in x sin <p 
in T sin 6 cos 6 — in, sin < 



(45) 



We define 7 = cos + in z sin </> and j3 = n x sin 0, and we 
write 



7 i/3 
e~ iGf 
iP 7* 



(46) 



Note that n x and n z are real-valued, hence P is real- 
valued, however 7 is complex. Evaluation of Eq. (38) 
yields 

Vi(t) = 

(7a + i/3b)e lGe 

( 7 *a-^)e~ iG£ (76 - i/3a)e~ lGe 

( 7 *6 + i/3a)e iG£ 

(47) 



The complex matrix elements of (47) are usefully ex- 
pressed as 



Then we have 



u = (ja + i/3b)e iGe , 
v= (jb-if3a)e- iGe . 



Vi(l) 




(48) 



(49) 



where |u| 2 + \v\ 2 = a 2 + b 2 = (a 1 /A£ v ) 2 = rj 2 . 

Heuristic construction for Ui(t, 0) 



D. 



Exact solution of wave equation (37) requires a time 
evolution operator Ui(t, 0) which solves Eq. (40) for the 
case wherein Vi(t) is a function of "live" variable t. A 
formal solution is provided in principle by the Dyson se- 
ries. In practice, the series is always truncated at low 
order; Sec. Ill of Ref. [5T] gives a clear discussion. As an 
alternative approximation which retains the perturbative 
expansion structure of the Dyson series, we introduce the 
exponentiation of Vj (£) as a heuristic form: 



U\ est \t = l,0>) 



-iVi(i)t 



(50) 



To obtain the matrix representation of Uj CSt ^ (£, 0) we 
take a brute force approach, rather than e.g., harness- 
ing the Cayley-Hamilton theorem [T^]. We observe that 



(v,y 



\u\ uv 

\u\ 2 + \v\ 2 



(51) 



(uvY 







Furthermore, \ VjJ 
by induction (Vi^j 



= V 

n— odd 



(Vi), (Vi) =v 2 (vi) 

= rf^Vi, and (V/) 



and 



n— even 



TJ 



Vi) . Thus 



-Wit 



E 

n=0 



(-iVii) n 




(52) 



(1 — cos r)£) 



.Vi . t 
1 — sm rj£. 

V 



The time evolution operator for neutrino propagation in 
our Interaction picture obeys a matrix identity reminis- 



cent of that for rotations generated by Jy 

Equation (52) yields an explicit representation for 
e -i%l, We write e = r/£, ? 



use (1 — cos 9 
of Eq. pj} is 



u/v, 



v/rj, and 



2 • sin I . Then the evolution operator 



l>j es4) (£,0) 



l-2|u| 2 sin 2 § -iusin 



9 — 2wvsin 



(53) 



-2(uv)* sin 



2 9 
2 



cos tt 
-iv* sin( 



1 



2|U| 2 sin 2 f 



We proceed to work our way back, first to the three- 
neutrino propagation basis and then to the neutrino fla- 
vor basis, wherein the matrix elements of U^(£,0) cor- 
respond to the possible neutrino oscillation amplitudes 
for 3-flavor mixing. The return of the time evolution op- 
erator to propagation basis requires that we calculate 



U (v) (£,{)) = 



iH, 



:) ' ■U ( I est \l,Q). 



(54) 



The required matrix forms are (53) and the adjoint of 



( 46 1 . As a preliminary to this multiplication, we define 



some compact forms. For the diagonal elements we define 

(55) 



D„ = 1 - 2|w| 2 sin 2 -, d = cos0, 



D„ 



l-2|?7rsin 2 



For the off diagonal elements, we define 



-2uv sin 



(56) 



Then we write 

jj(est) 



(£,0) = 





—ie p 


-ie* 


d —ik 


P* 


-ik* D v 



(57) 



Proceeding with the evaluation of (54): 

U {p \£,Q) = 

("f*D u -i/3p*) (7*(-ic) - pk*) (j*p-ipD v ) 
{-ie*) e lGl d e lGe (-ik) e lGl 

( 7 P* - ipD u ) (7(-***) - Pe) (7 A, - ijSp) 



Finally, we return to neutrino flavor basis via 

u^\e,o) = Ai(023) • h CP ■ u [p) (e,o) ■ L Scp ■ R[(e 23 ). 

All of the physically relevant neutrino flavor oscillation 
amplitudes are contained in the matrix U^(£,0); how- 
ever our focus here is upon the element u[%* ■ 

E. — > v e oscillation amplitude 

Element provides the v e appearance amplitude 

from an initial beam of neutrinos: 

U[f = A(y M -► v e ) = onUjg + s 23 ug>e- i5 ° p . (58) 
Unwinding notations and rearranging, we obtain 



Now ab/rj 1 — — C13S13 and a 2 /rf — c 2 3 . Recall that 
7 = cos 4> + in z sin <f>, and that (3 = n x sin cf> = if sin <fi = 



sin26» 13 • fffi. Working within the bracket of Eq. (|63|), 



we separate the real and imaginary pieces: 



T 3 =2s 23 



sin 26* 



13 



cos + i 



sin 2^13 2 
- n z + c^ 3 n x 



Then 



T 3 = sin 26»i3 s 2 3 



cos{N£) + iF A 



sin(iW) 
AL,N 



I"., I 



which is identical to Eq. (JsJ) . 



^ -> i/ e ) = (-i)«28^e- Woi ' 

+ ( — i)c23 [7*m — i/3CS*] • sin0 

Q 

+ 2s 23 [if3\v\ 2 - j*uv] ■ sin 2 - • e~ i5ci 



(59) 



We identify the three terms of Eq. ( 59 1 using the ex- 
pression Atyn -> v e ) = Ti + T 2 +T3. For Ti we have 



Ti = (-i)s 23 /3e lfcp = (-i)s23n a sm^- e 



(-i) sin 20i3 S23 



sin 



(60) 



which coincides with Eq. 
For T 2 we have 

T2 = — 1023 [7*^ — sin0 

Assembling the various factors, 

T 2 = -i— [7* (70 + 0b) - i0y*b + P 2 a] 
i] J 

= -ic 23 • - • sin(0) • e lGt , 
V 

which is the same as Eq. 

For the remaining term T 3 we have 



(61) 



sin f ■ e 



iGl 



T 3 = 2s 2 3 [*/?M 2 - 7*™] • sin" - ■ < 



2 v _ p -iScf 

2 



(62) 



V. CONCLUSION 

Having presented and subsequently derived each of the 
probability terms individually, we conclude by stating the 
entire six-term formula for Viy^ v e ) : 



, ■ x2 2 sin 2 (£>A) 
(sin20i 3 ) 2 s\ 3 - y D2 ' 

+ sin 20i3 c i3 sin 2023 ' sin(a' A) 



sin(DA) 
D 



[cos A' • cos 8c p — sin A' • sin 6c p] 

c 2 13 c 2 23 -sm 2 (a'A) 

2 „ . o /VAN sin 2 (L>A) 
2 sin 20i3 sin 20i 3 s 2 3 F A ■ sin 2 ' ] 



sin 20i3 c i3 sin 20 2 3 • sin(a' A) • sin 



2 / D 2 
a' A' 



{cos(£>A) • sin(A' + S CP ) ~ F A Sm{ ^ A) . CO s(A' + S C p)\ 



-sin 2 20i3 & • sin 4 f'^) . jcos^A) + F 2 ■ 

(64) 



As discussed in Sec. HID the three leading terms of Eq. 

641) are reminiscent of the perturbative expansion of Eq. 

TIT These three terms account for essentially all of the 
Vft — > v e appearance probability. 



Reduction within the bracket leads to 

,2 " 



T 3 = -2s 



2:i 



ab 
i] 2 



ij3- 



sin 2 - • e- lScr 
2 



(63) 
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